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Abstract 

Transformation elasticity, by analogy with transformation acoustics and optics, con- 
verts material domains without altering wave properties, thereby enabling cloaking and 
related effects. By noting the similarity between transformation elasticity and the theory 
of incremental motion superimposed on finite pre-strain it is shown that the constitutive 
parameters of transformation elasticity correspond to the density and moduli of small- 
on-large theory. The formal equivalence indicates that transformation elasticity can be 
achieved by selecting a particular finite (hyperelastic) strain energy function, which for 
isotropic elasticity is semilinear strain energy. The associated clastic transformation is 
restricted by the requirement of statically equilibrated pre-stress. This constraint can 
be cast as trF = constant, where F is the deformation gradient, subject to symmetry 
constraints, and its consequences are explored both analytically and through numerical 
examples of cloaking of anti-plane and in-planc wave motion. 

1 Introduction 

The principle underlying cloaking of electromagnetic and acoustic waves is the transformation 
or change-of- variables method [H [2] whereby the material properties of the cloak are defined 
by a spatial transformation. While the first applications were to electromagnetism, e.g. [3], it 
was quickly realized that the same mathematical methods work equally well in acoustics [U [U 
[6]. The fundamental identity underlying electromagnetic and acoustic transformation is the 
observation that the Laplacian in the original coordinates maps to a differential operator in the 
physical coordinates that involves a tensor which can be interpreted as the new, transformed, 
material properties [7] . The equivalence between the Laplacian in the original coordinates and 
the new operator involves an arbitrary divergence free tensor [8] , implying for the acoustic case 
that the transformed material properties are not unique. For a given transformation function 
one can achieve acoustic cloaking by a variety of materials, ranging from fluids with anisotropic 
inertia, to quasi-solids with isotropic density but anisotropic stiffness [HI [9] . Non-uniqueness 
of the material properties does not apply in the electromagnetic case, where, for instance, 
the permittivity and permeability tensors must be proportionate for a transformation of the 
vacuum. 

A crucial aspect of transformation optics and acoustics is that the mapped properties 
correspond to exotic material properties far removed from the realm of the original material. 
This aspect is accentuated in transformation elasticity. In the first study of transformation 
theory to elastodynamics, Milton et al. [TO] concluded that the transformed materials are 
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described by the Wihis model. This constitutive theory for material response is dispersive, 
involving coupling between stress and velocity, in addition to anisotropic inertia |lHll2j. Brun 
et al. [13] considered the transformation of isotropic elasticity in cylindrical coordinates and 
found transformed material properties with isotropic inertia and elastic behavior of Cosserat 
type. The governing equations for Cosserat elastic materials |14j are the same as those of 
"normal" linear elasticity except that the elastic moduli do not satisfy the minor symmetry, 
i.e. Cjfj^i / ^ijki (although they do satisfy the major symmetry C^f- = CfJ|;). This implies 
a non-symmetric stress, cr^ ^ cr which depends not only on the strain e (the symmetric part 
of the displacement gradient) but also on the local rotation |(Vu — (Vu)*). 

A thorough analysis of transformation theory for elasticity [15] indicates that, as in acous- 
tics, the range of mapped material properties is highly non-unique, thus explaining the diver- 
gence in the previously obtained results [inillS]. The transformed elastodynamic constitutive 
parameters may be characterized through their dependence on (i) the transformation (map- 
ping function) and (ii) on the relation between the displacement fields in the two descriptions, 
represented by matrices: F, the deformation matrix, and A, respectively. It was shown [15] 
that requiring stress to be symmetric implies A = F and that the material must be of Willis 
form, as Milton et al. |10j found. Setting A = I, on the other hand, results in Cosserat 
materials with non-symmetric stress but isotropic density, as found by Brun et al. [13] and 
by Vasquez et al. [16]. An alternative approach to transformation elasticity has been pro- 
posed that employs inextensible fibers embedded in an elastic material \n\ [T8] . This has the 
advantage that the effective material has isotropic density and retains both the minor and 
major symmetries of the stiffness tensor. Despite the non-unique nature of transformation 
elasticity the materials required are, in all cases, outside the usual realms of possibility. 

In this paper we consider a class of materials displaying non-symmetric stress of the type 
necessary to achieve elastodynamic cloaking. Effective moduli with the major symmetry 
{Cijki = Ckiij) that do not display the minor symmetry {Cijki / Cjiki) are found in the 
theory of incremental motion superimposed on finite deformation [19]. We take advantage of 
the similarities between transformation elasticity and small-on-large motion in the presence of 
finite pre-strain. The starting point is the formal equivalence of the constitutive parameters 
of transformation elasticity with the density and moduli for incremental motion after finite 
prestress. This offers the possibility for achieving elasticity of the desired form by proper 
selection of the finite (hyperelastic) strain energy function. Such an approach has been shown 
to be successful in the context of antiplane (SH) wave motion. By using the neo-Hookean 
strain energy for incompressible solids and applying a radially symmetric cylindrical pres- 
strain, Parnell [20] and Parnell et al. [21] showed that the resulting small-on-large equations 
are identically those required for cloaking of the SH wave motion. Here we consider the 
more general elastic transformation problem, including but not limited to SH motion. We 
show that the form of the finite strain energy is restricted in form for isotropic elasticity. 
The equivalence between the transformation and the finite pre-strain also limits the type 
of transformation possible. This contrasts with the acoustic and electromagnetic problems 
for which the transformation is arbitrary. The elastic transformation is restricted in form 
because the pre-stress must be statically equilibrated, implying that the transformation must 
satisfy a partial differential equation. We show that this constraint can be cast as trF = 
constant, (subject to symmetry constraints), and explore its consequences both analytically 
and numerically. 

We begin in ^ with a review of transformation elasticity and of incremental motion 
superimposed on finite pre-strain, emphasizing equivalence of the theories. The form of the 
finite strain energy necessary to achieve transformation elasticity is deduced in ^ and in 21 
the constraint on the deformation for isotropic elasticity is derived. ^ presents a detailed 
example of the type of radially symmetric finite pre-strain possible for isotropic elasticity. 
These analytical results are extended and illustrated in ^ through numerical examples of 
cloaking of anti-plane and in-plane wave motion. Conclusions are presented in ^ 
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2 Background equations 



We first review the theory of transformation elasticity for linearly elastic solids, and then 
consider the separate theory for incremental deformation in finite elasticity. 



2.1 Review of transformation elasticity 
2.1.1 Transformation notation 

A transformation from the virtual configuration, f^o, to the present configuration ft (also 
known as the physical or current domain) is described by the mapping from ^ € JIq to 
X G ri. Component subscripts in upper and lower case (/, J, . . . , i, j, . . .) are used to distin- 
guish between explicit dependence upon ^ and x, and the summation convention on repeated 
subscripts is assumed. The transformation or mapping is assumed to be one-to-one and in- 
vertible. Perfect cloaking requires that the transformation is one-to-many at a single point 
in Qq. This can be avoided by considering near-cloaks, where, for instance, a small hole in 
r^o is mapped to a larger hole in Q. 

The transformation gradient is defined as F^'^^ = V^x with inverse F^'^^ = V^, or in 

component form pj^j'' = dxi/d^i, pj^j'' = d^i/dxi. The Jacobian of the transformation is 
Jo = detF^'^-'. The infinitesimal displacement u^'^)(^,i) and stress a^^\$,,t) are assumed to 
satisfy the equations of linear elasticity in the virtual domain: 

div5<T(0) =pou(o), aW = C(0)V5u(0) infio, (1) 

where pQ is the (scalar) mass density and the elements of the elastic stiffness tensor satisfy the 
full symmetries c}j^^ = Cj^jj^j^, cfjj^^ = C^\jj; the first identity expresses the symmetry 
of the stress and the second is the consequence of an assumed strain energy density function. 

Particle displacement in the transformed domain, u(x, t), is assumed to be related to the 
displacement in the virtual domain by the non-singular matrix A as 

u(o)=A*u {uf^=Auu,). (2) 

The choice of the transpose. A* in eq. ([2]), means that the "gauge" A and the transformation 
gradient F^*^^ are similar objects, although at this stage they are unrelated. Neither A or F^*^^ 
are second order tensors because of the fact that they each have one "leg" in both domains. 
Milton et al. [10] specify A = F^*^^ since this is the only choice that guarantees a symmetric 
stress (see §2.1.2p . Identifying [10] d^ and dx with u^^^ and u, respectively, and using 

d^ = F*-") dx would lead one to expect A* = F^''^ . However, the displacements are not 
associated with the coordinate transformation and F^*^^ and A are independent quantities. 



2.1.2 The transformed equations of elasticity 

Under the transformation (or change of coordinates) ^ — )• x the equilibrium and constitutive 
relations ([T]) transform to jl5] 

crij,i = Pj, (Tij = Cfj^iUi^k + Sfji ill, Pi = Siji Uj^i + p'ji iij, in fio, (3) 

with parameters C^, S'^^ and p'^^ defined as follows in the Fourier time domain (dependence 
g-ri^t understood but omitted) 

Stfi = {-iu:)-^JoCf]KLQi3ijQkiKL,k, (4) 

Pjf = PoJo'^^jKAiK + {-i^^)~^JoC\°]j^^Qijij^iQklKL,k, 
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where Qijij = Jq ^F^j^ Ajj. The elastic moduh and the density satisfy the symmetries 

^eff _ /^cS off _ eff ft:\ 

'-'ijkl — ^kliji Pjl — Plj > l-^J 

although these are not the full symmetries for the Willis constitutive model (which requires 
the additional "minor" symmetry C^^^^ = Cjfj^i). Equations ©-([I]) are the fundamental result 
of elastic transformation theory |15j . 

The absence of the minor symmetries under the interchange of i and j in C^^^i and 
of dH) implies that the stress is generally asymmetric. Symmetric stress is guaranteed if 
QijiJ = QjiiJ, and occurs if the gauge matrix is of the form A = (^F^*^), for any scalar C / 0, 
which may be set to unity with no loss in generality. This A recovers the results of Milton 
et al. [TO] that the transformed material is of the Willis form. As noted in [10], it is the only 
choice for A that yields symmetric stress. 

The equations in the transformed domain, which is the physical realm, clearly display a 
great deal of non-uniqueness, corresponding to a vast realm of possible material properties. 
Our preference is for non-dispersive (i.e. independent of frequency) materials, in particular, 
those with the least "unusual" properties, so that they can conceivably be related to actual 
materials. In this regard, isotropic density is achieved by taking the constant matrix A 
proportional to the identity, A = (^I, with C = 1 without loss of generality. In this case 
p'^ff = S^^ = 0, with non-dispersive density and elastic moduli given by 

^cff _ T-l ^ ^cfT _ T-l p(0) p(0) ^(0) ff.^ 
P - Jq PO, ^ijkl--Jo ^il ^kK^IjKl- W 

The equations of motion in the current domain are then 

{C^^,iUi,k),i = P'''uj. (7) 

The effective moduli of satisfy the major symmetry ([5|)i but Cj?||; 7^ ^jfkv indicating 
a non-symmetric stress. Departure from symmetric stress is possible in continuum theories 
such as Cosserat elasticity and micropolar theories of elasticity. Another context admitting 
non-symmetric stress is the theory of small-on-large motion, described next. 



2.2 Small-on-large theory 

The solid material is considered in two distinct states: first, the reference configuration of 
the solid under zero strain, Qi, and secondly the current state of the material, which is 
again identified with Q. The hyperelastic theory of small motion superimposed upon large 
depends upon the initial finite, i.e. large, static pre-strain which maps X G to x G ri. The 
subsequent small motion is defined by the dynamic mapping X — t- x -|- u(x, t). The following 
theory assumes u and the associated strain are sufficiently small that tangent moduli can be 
employed to derive the linear equations of motion for the small-on-large motion [19] . 

Towards that end we introduce the deformation gradient of the pre-strain F = VjfX with 
inverse F~^ = VX, and Jacobian J = detF. The polar decomposition is F = RU = VR, 
where R is proper orthogonal (RR* = R*R = I, det R = 1) and the tensors U, V G Sym"'" 
are the positive definite solutions of = C = F*F, = B = FF*. The material is assumed 
to be hyperelastic, implying the existence of a strain energy function W per unit volume from 
which the static Cauchy pre-stress is defined as 

afp = J-^Fi^dW/dFj^. (8) 

The assumed dependence of W on the deformation F, along with the freedom to change the 
current coordinate basis (which has nothing to do with transformation!) implies that W must 
depend upon QF for any orthogonal Q, and taking Q = R* implies the dependence W = 
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yy(U). Assuming the density in the reference configuration is pr, the governing equations for 
subsequent small-on-large motion u(x, t) then follow from the well known theory |19j as 



{•Aoijkiui,k),i = puj,tt, (9) 



where 



_i _i d'^yv 

p=J Pr, Aoijkl = J FiaFkp AajpU Aaj/^l = — 777;— (= (10) 

3 Potential strain energy functions 

Our objective is to find possible hyperelastic solids, i.e. strain energy functions W such that 
the equations for small-on-large motion are equivalent to those required after transformation 
of a material assumed to be homogeneous with properties {po, Cj^^j^j}. 

The connection between the transformation and the small-on-large theories is made by 
first identifying the displacement fields as equivalent, u(x, t) = u(x, t), and then requiring 
that the equations of motion ([7]) and ^ are the same. The latter is satisfied if 

p=7P''^ Ao^,kl=7c:^kl, (11) 

for some positive constant 7. Hence, 

J- Vr = 7^0"' Po, F,o.F,p Aaj^i = 7 Jo' F^^^FlfjI. c\%. (12) 

The reference density pr can then be chosen so that 7 = 1, and eq. (|12p then implies that the 
hyperelastic material is defined by 

Pr = PoJq'J, Aaj(3l = Jq'J F~l F^f F^^f'^^ ^fjKV (1^) 

Equation (|13p 1 is automatically satisfied if the transformation and the finite deformation are 
related in the following manner: 

F= (Wpo)'^'fWG-\ 5 = detG, (14) 



for some non-singular G. Equation (I13p 9 combined with the expression for Aajpi in eq. (flO 
yields a second order differential equation for the strain energy function. 



dFj^dFip 



Recall that pq and C'/jXi constant, but at this stage the remaining quantities in (|15p . 
i.e. Pr and G, are not so constrained. The density in the reference configuration could be 
inhomogeneous, pr = pr(X.)- In that case (jlSp would not have a general solution for W unless 
G also depends upon X in such a manner that the right hand side is independent of X. 
This suggests that permitting pr to be inhomogeneous does not provide any simplification, 
and we therefore take the reference density to be constant, although not necessarily the same 
as Pq. The quantity G could, in principle, be a matrix function of F, but this makes the 
integration of (jl5p difficult if not impossible. We therefore restrict attention to constant G. 
Consideration of the important case of isotropic elasticity in fg] indicates that the degrees of 
freedom embodied in G do not provide any significant additional properties, and therefore 
for the remainder of the paper we take G = I, and set pr = po with no loss in generality. In 
this case the solution of (jlSp such that W = under zero deformation (F = I) is 

W = liF,^ - 5,^){Fip - 6ip) C%,. (16) 
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Equation (|16p provides a formal solution for W, one that is consistent with (|15p . However, 
the dependence of W in (|16p upon F points to a fundamental difficulty, since the strain energy 
should be a function of U. The two are not equal in general, unless 

R = I ^ F = U = V. (17) 

We henceforth assume (jl7p to be the case: that is, we restrict consideration to deformations 
that are everywhere rotation-free. Equation then suggests the following possible form of 
the finite strain energy 

W = ^EjaEi^C^^^^i where E = U - I. (18) 

Although this has realistic dependence on U, it will not in general satisfy eq. (jlSp . i.e. 
d'^yV / dFjadFip 7^ ^ajpv return to this crucial point for isotropic elasticity in ^4.31 where 
we demonstrate that eq. (fT5]) is satisfied by the isotropic form of ([18]) under additional con- 
ditions. Note that the strain measure E, which is sometimes called the extension tensor., has 
as conjugate stress measure Sq = 9W/5E = ^(SU + US) where S = JF~-'^<t(F*)~^ is the 
second Piola-Kirchhoff tensor \12\ §2.5]. 

We restrict attention henceforth to the case of hyperelastic materials that are isotropic in 
the undeformed state. 



4 Isotropic elasticity 

4.1 Semilinear strain energy function 

The initial moduli are C^^^pi ~ ^o^aj^pi + Mo('^a/3<Jj7 + ^la^jp) with original Lame moduli 
^0 > 0, Aq and Poisson's ratio v = Ao/[2(Ao + /^o)] £ We consider the isotropic 

version of the hyperelastic strain energy in (jl8p . 

W = y (tr E)2 + ^0 tr(E)2 = ^(n - 3)^ + /.o((n - 1)' - 2(^2 - 1)), (19) 

with the latter expression in terms of two of the three invariants of U: ii = Ai + A2 + A3, 
i2 = A1A2 + A2A3 + A3A1 where Ai, A2, A3 are the principal stretches of U. Materials with 
strain energy (jl9p have been called semilinear [23] because of its relative simplicity and the 
linear form of the Piola-Kirchhoff stress T^, related to the Cauchy stress by <tP''° = J~^FT^; 
thus Tr = 2^qF + (Ao(trE) — 2//o)R- John [H] proposed the strain energy ([T9]) based on the 
explicit form of its complementary energy density in terms of T^j, a property also noted by 
others |25| [26] . The semilinear strain energy is a special case of the more general harmonic 
strain energy function [23]. Plane strain solutions for harmonic strain energy are reviewed in 
[271 §5.2]. Sensenig |28] examined the stability of circular tubes under internal pressure, while 
Jafari et al. [29] considered both internal and external pressure loading. The latter study has 
implications for the stability of the pre-strain developed here, see ^5.2.31 
The pre-stress follows from ([8|) as 

= J-i [Ao(n - 3)V + 2^o(v2 - V)] . (20) 

It is emphasized that we are restricting attention to deformations with U = F = F*, so 
that the Piola-Kirchhoff stress is also symmetric with T/j = Ao(tr E)I + 2^qEi. Applying the 
equilibrium equation for the finite deformation, 

DivTfl = XoXa^aj + '^^J■OXj,aa = 0. (21) 
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We seek solutions with symmetric deformation gradient, Xj^a = ^a,j, and hence Xa^aj = Xj,aa- 
Consequently eq. (j2ip is satisfied by finite deformations satisfying either of the equivalent 
conditions 

Xj,aa — and Xa^ctj — 

0. Thus: 

Since the two partial differential equations in (j22p are the same we need only seek solutions 
of one. Focusing on x^^aj = we conclude that the most general type of deformation x(X) 
is described by 

Divx = c (= constant > 0), where Vxx = (Vjsfx)"^. (23) 

The appearance of the positive constant of integration in (j23P i means that the sum of the 
principal stretches is fixed, 

Ai + A2 + A3 = c {ii=c). (24) 

Further implications of the general solution (j23p for a material that is isotropic in its un- 
deformed state is explored in greater detail in the next section. For now we note that the 
pre-stress follows from ([20]) and (p^ as 

= 2/.0 (V^ - V + ^^^V) . (25) 



4.2 The limit of z/ = | 

It is of interest to consider the limit of the isotropic solution for = |. By assumption the 
prestress must remain finite. Consequently, using eq. (j25p . it becomes clear that in the limit 
as — >• ^ the constant of integration c = 3, i.e. 

W = /iotr(U2-I), <tP'-^=pJ-^V + 2^oJ-'(V2- V),] ^ 1 . ^ 

_ > for = - (26) 

Divx = 3, F = F^(=U = V), J 2 

where the scalar p(X) defines the constraint reaction stress (the factor is included for 
later simplification). The latter arises from the limiting process of z^ — )• ^ in eq. (j25p . and 
has also been shown to be the unique form of the reaction stress for the constraint tr V = 3 
|30| . Note that in writing <tP''° in (j26p we maintain a term proportional to V in the second 
term rather than incorporating it with the constraint term. This form is consistent with the 
requirement that p = and hence a^'^'^ = in the undeformed state x = X. The equilibrium 
equation for the pre-strain follows from eq. (j26p as V xP + 2/ioV^x = 0, and since Vj^x = 
(see eq. (|22p) it follows that p = constant. 

Several aspects of (j26p are noteworthy. The limit of z/ = ^ is usually associated with 
incompressibility, i.e. the constraint J = 1 or equivalently i^, = A1A2A3 = 1, although the 
reason underlying this identification originates in linear elasticity and is therefore by no 
means required. Strictly speaking, the isochoric constraint is = 1 conserves volume under the 
deformation. Here we find that z^ = ^ implies the kinematic constraint on the deformation 
that ii = Ai + A2 + A3 = 3. The latter is associated with the notion of incompressibility 
in linear elasticity in the form ti E = Q and in the present context can be viewed as a 
"semilinear" feature, in keeping with the descriptor [23] for the strain energy function ()19p . 
The kinematic condition, Div x = 3 or equivalently 

Ai + A2 + A3 = 3 (trV = 3), (27) 

has been called the Bell constraint [30] by virtue of the fact that Bell [31] [32] showed it to 
be consistent with numerous sets of data for metals in finite strain. Solids satisfying this 
constraint have been called Bell materials [30]. In contrast to the constraint A1A2A3 = 1 
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it can be shown that volume decreases for every deformation of a Bell material, and hence 
isochoric deformations are not possible j30j . 

Another feature of the i' = ^ limit is that the strain energy in (j26p has the functional 
dependence W = /xo(^i + ^2 + -^i ~ 3)- interesting to compare this with the strain energy 
for a neo-Hookean solid, Wnh = ^(-^1 + ^2 ~^ ^3 ~ ^ associated with incompressibility (i.e. 

= 1). Both strain energies reduce to the incompressible form for linear elasticity, and the 
factor i appearing in W but not in Wnh can be attributed to the different constraints in 
each case {ii = 3 or 13 = 1). Parnell [2U] and Parnell et al. [21] considered neo-Hookean 
materials in the context of transformation elasticity for isochoric deformation. The present 
results indicates that the same form of the strain energy but with a different constraint yields 
a quite distinct class of volume decreasing deformations. This aspect will be examined further 
in the next section in terms of specific examples. 

4.3 Consistency condition 

It remains to show that the semilinear strain energy ()19p satisfies 

■^ajl3l = X^SajSjil + Hoi^al^Sjl + ^la^jis), where Aajl3l = ^ ^ • (28) 

Since the moduli .A. are isotropic, it it sufficient to show the equivalence in any orthogonal 
system of coordinates. We choose the principal coordinate system, in which the non-zero 
components of A, for isotropic elasticity satisfy [191 eqs. (3.31)-(3.34)] 







= y^^v 




(29a) 


Aijij 


— Aijji 


_Wi + Wj 
Xi + Xj 




(29b) 


Aijij 


+ Aijji 


_Wi- Wj 
Xi — Xj 

= y^ii - Wij, 


« / h \ / Aj, 


(29c) 


Aijij 




^ / h \ = Aj, 


(29d) 



where Wj = dW/dXi, Wij = d'^W/dXidXj , i,j G {1,2,3} with no summation on repeated 
indices. Using W and c as defined in eqs. (fT9l) and (i24|l gives 

= Ao (c-3) + 2;uo {Xi-l), 
= Ao + 2^0, = Ao, j. 
These satisfy (f29iD . ([29c]) and (f29d]l . The remaining conditions (I29bp become 

Wi + W,=0 ^ (Ao + ^o)(c-3)-Mo(Afc-l) = 0, i^j^k^i. (31) 

Equation (|3ip constitutes three conditions, which taken together imply the unique but trivial 
solution Aj = 1, i G {1, 2, 3}, i.e. zero pre-strain. We avoid this by restricting attention to two 
dimensional dynamic solutions only, either in-plane (P/SV) or out-of-plane (SH) motion. 

4.3.1 In-plane (P/SV) motion 

The small-on-large displacements for in-plane motion are of the form ui{xi,X2,t), U2{xi,X2,t), 
U3 = 0. The condition ()3ip then only needs to be satisfied in the single instance i,j = 1,2, 
implying that the out-of-plane extension is related to the sum of the in-plane extensions by 

A3 = l-^(Ai + A2-2). (32) 

Since A3 is strictly positive, this places an upper limit on the sum of the in-plane extensions: 
Ai + A2<2(l + i/). 
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(30) 



4.3.2 Out-of-plane (SH) motion 



The out-of-plane SH motion is of the form ui = U2 = 0, U3{xi,X2,t). The requirement now 
is that ^1313 and ^2323 are both equal to ^0 i^i order to recover the out-of-plane equation of 
motion and associated tractions. Using ()29bp and ()29cp 

-41313 - MO = (^^) [c - 3 - (1 - 2u){X2 - 1)] , 

7 + ^=^ (33) 

^2323 - M = (f^) [C - 3 - (1 - 2z.)(Ai - 1)] , 

M + A3 

where c is the constant from eq. ()24p . In this form it is clear that if 7^ ^ then in-plane 
pre-stretches must be the same, Ai = A2 = 1 -|- (c — 3)/(l — 21/) , and therefore all the stretches 
are constant (since c is a constant). This rules out the possibility of SH cloaking since we 
require that the in-plane pre-strain be inhomogeneous. However, if both ^ = \ and c = 3 
simultaneously hold, then .A1313 = ^2323 = /^o for inhomogeneous and unequal in-plane 
stretches Ai and A2. We are therefore led to the conclusion that SH cloaking requires a 
separate limit of the semilinear strain energy, one satisfying the Bell constraint ()27p for which 
the strain energy and stress are given by (|26p. Note that we do not get the neo-Hookean 
strain energy in this limit. 



5 Applications to isotropic elasticity 
5.1 Radially symmetric cylindrical deformations 

Consider deformations that are radially symmetric, r = r{R), 9 = Q, in cylindrical coordi- 
nates (r, 9, X3) and {R, G, X3). The stretch in the xs-direction is assumed fixed, A3 = constant. 
The deformation gradient for r = r{R) is irrotational with 

(F* = ) F = Xrlr + Xele + hh, K = r\ >^e = (34) 

where 1^ = (8> , I^i = eg eq and I3 = 63 Cgi 63. The condition (|24p implies that the sum 
of the in- plane principal stretches is constant, say cq, and the constraint ()32p relates this to 
c of eq. dlZl), 

(1 - 2v)co + 2vc = 2(1 + v) where cq = A^ + Ag. (35) 
Equation (|23p for x reduces to an ordinary differential equation for r{R), 

r' + Jl = CO, (36) 

with general solution 

r = -^R + ciR~^, ci = constant. (37) 



Note that the free parameter cq may be expressed in terms of either c or A3, using eqs. (j32p 
and ([35|) . Using eqs. ([25]) . ([36]) and ([Mj) it follows that the principal stretches and stresses 
for the radially symmetric cylindrical configuration are 

A, = 2- Ae + 21/(1 -A3), Ae = ^, A3, 

(38) 

ar = T^(A.-A,), ar = T^(A, - A.), a^f = |^(1 + z.)(A3 - 1). 

A3A0 A^Ar ArAg 

Note that 

^ ^ 1, ^ as r ^ 00 iff A3 = 1 {^cq = 2). (39) 

R 
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5.2 Two types of cloaking 
5.2.1 Conventional cloaking (CC) 

The conventional concept of a cloaking material is that it occupies a finite region, in this case, 
the shell R £ [A, B] which maps to an equivalent shell in physical space with the same outer 
surface and an inner surface of radius larger than the original, i.e. r G [a, -B], a € {A,B). 
Applying (j37p with the two constraints r{A) = a and r{B) = B yields 



R+{a-A) 



B^2 



(I) 



1 



R 

A' 



Re[A,B], 



(CC) 



(40) 



which specifies the previously free parameter A3 (also c and cq) as 

1 



A,3 



1-i 

V 



a_ 

A 



1 



< 1. 



(41) 



The constraint A3 > therefore sets a lower limit in the permissible value of the outer radius 
as 



B> A 1 + 



1 , a 
u^A 



1 



1/2 



(42) 



The mapping (j40p must also be one-to-one within the shell with A,. = r' > 0. This means 
that there should be no zero of r' = for R £ [A,B]. The convex nature of the solution 
(jlO|) implies there is only one zero, say at R = Rq. Since sgnr' =sgn(i? — Rq), it follows that 
Rq < A must hold. Noting from (p7|) that Rq = 2ci/cq, and using (jJO]) to infer cq and ci, the 



condition Rq < A? becomes 



r A 2 

a < amax = 2^/^l + (— ) 



(43) 



The magnification factor > 1, which measures the ratio of the radius of the mapped hole 
to the radius of the original one, is therefore bounded according to 



A 



< 2. 



(44) 



In order to achieve a reasonable degree of cloaking one expects that the magnification factor is 
large, so that the mapped hole corresponds to an original hole of small radius and hence small 
scattering cross-section. The limitation expressed by (|44|) therefore places a severe restriction 
on the use of the hyperelastic material as a conventional cloak. Note that the upper limit on 
a in (fl3l) is not strictly achievable because a = Omax implies Xr{A) = r'{A) = and hence the 
principal stresses cr^^'^, cr^z become infinite at r = a. 

The hyperelastic mapped solid has other aspects that further diminish its attractiveness 
as a conventional cloaking material. Specifically, a non-zero traction must be imposed on 
both the outer and inner boundaries to maintain the state of prestress. Noting that the 
radial stress is 

2/10 R \ 'T 



a. 



pre 



{R) 



A3 



R-^ 



v(\ - A3) 



with r given in (j40p . yields 



ar(A) = ar(i?)-^(i--), 

A3 a 



CTr(i?) = f^Kl-A3). 



(45) 



(46) 



The necessity of the inner traction at r = a is a reasonable condition, but the requirement 
for an equilibrating traction at r = i? is physically difficult. One way to avoid this is to let 
— 7- 00, considered below. 
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Examples of the radial deformation are given in Figure [TJ 



r 



0.2 



0.8 
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0, 
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0.6 



0.8 



R 



Figure 1: The transformed radius r as a function of R for 2D radially symmetric pre-strain, 
from eq. (j40p . The seven curves correspond to ^4 = 0.1, 0.2, . . . , 0.7 with i? = 1 in each case. 
The value of the mapped inner radius for each curve is a = r{A) = O.Oomax where Omax 
is defined in eq. ()43p . The dashed line indicates r = R. Mappings that lie above this line 
represent spatial compression, applicable to cloaking. 



5.2.2 Hyperelastic cloaking (HC) 

The hyperelastic material is now considered as infinite in extent. The cloaking effect is caused 
by allowing a radially symmetric hole in the unstressed configuration to be expanded under 
the action of an internal pressure to become a larger hole. We therefore require that the 
traction at infinity is zero, and that r/R tends to unity, so that (|39p applies. Then setting 
the mapped hole radius to r{A) = a {> A) implies the unique mapping 



This deformation is simply the limiting case of (j40p for i? — )■ oo. Note that the restriction (j44p 
still applies to the magnification factor in order to ensure A,. > for r > a. The traction 
at the inner surface is a pressure which follows from (j46p in the limit — )■ oo, A3 — )• 1, as 



It is interesting to note that the internal pressure is independent of the Poisson's ratio and 
it is therefore the same as pi^ found by |20j considering SH incremental motion. 

5.2.3 Stability of the pre-strain 

Jafari et al. [29] examined the stability of a finite thickness tube composed of material with 
harmonic strain energy, which includes semilinear strain energy as a special case. They 
showed that radially symmetric two dimensional finite deformations are stable under interior 
pressurization with zero exterior pressure. This implies that the finite pre-strain HC is stable. 
The stability of the CC deformation (|40l) does not appear to have been considered, and remains 
an open question. However, the stability of the HC, corresponding to S — > 00, means there 
exists a minimum i?min for which CC stability is ensured for all B > i?min- 



r = R + {a- A){A/R), Re[A,oo). 



(HC) 



(47) 




(48) 
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5.3 The limiting case when = 2 

In this limit the constraint (j27p apphes and the pre-stress for the radially symmetric defor- 
mation follows from (j26p with constant "pressure" p (see ^4.2p as 

= ^{K - 70), = Y^i^e - 70), (49) 

where the value of the constant 70 = (2/LIo) depends on the specified boundary conditions, 
and \r = dr/ di?, \q = r/R, with r{R) given by eq. ([37|) for cq = 2, For instance, in the 
case of hyperelastic cloaking as defined in ^5.2.21 we find, noting the result (p9]) . that p = 0, 
yielding the same interior pressure p\a as eq. (jl8|) . 
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Figure 2: SH wave field. Left: Total (top) and scattered (bottom) fields corresponding to an 
undeformed cavity with scaled radius KgA = 2tt. Right: Total (top) and scattered (bottom) 
fields corresponding to a conventional cloak generated via pre-stress where the scaled deformed 
inner radius is Kga = 27r and initial inner cavity radius defined by a = I3A where /3 ~ 1.863. 

6 Numerical examples 

We illustrate the above theory in the two dimensional setting where we consider wave scat- 
tering from a cylindrical cavity with and without a cloak where the cloak is a conventional 



12 




Figure 3: Scaled scattering cross-section ^shKs (left) and percentage reduction in scattering 
cross-section 7% by using a hyperelastic cloak (right) (with a = j3A where j3 k, 1.863.), both 
plotted against scaled cavity radius Kga for the SH wave case. The cross-section is plotted 
without (solid) and with (dashed) a hyperelastic cloak. A significant reduction in scattering 
is achieved by using a hyperelastic cloak. 



cloak created via pre-stress. We shall show that partial cloaking is achieved, in the sense that 
scattering is significantly reduced by presence of a cloak. We are not able to achieve perfect 
cloaking since the cavity has to be of finite radius initially and furthermore, the hyperelastic 
theory above restricts the expansion to be at most twice the initial radius, i.e. a < 2A. We 
consider two cases: horizontally polarized shear (SH) waves and coupled compressional/in- 
plane shear (P/SV) waves. We take B/a = 2 which upon using (j43p gives an initial inner to 
outer cloak radius ratio B/A = l/{2- a/S) 3.732 and a/ A = f3 = 1/(2(2 - V3)) ^ 1.866. 

For the SH and P/SV wave examples considered below, we use the description in Appen- 
dices lA.ll and IA.2I regarding scattering from a cylindrical cavity in an undeformed medium 
due to an incident field generated by a line source. In both cases considered, we assume that 
the line source is located at a distance Rq from the centre of the cavity with Rq/B = 2, it is 
of unit amplitude C = 1 and in the P/SV case it generates purely compressional waves. 

6.1 SH wave propagation 

In this case the shear wavenumber Ks of the medium is defined by K'^ = w^/c^ = pquP' / jiQ 
where pQ is the density of the medium in the undeformed configuration. We use the solution 
in Appendix lA.ll to solve the corresponding (conventional, pre-stress) cloak problem, the 
difference arising merely due to the modified argument due to the hyperelastic deformation 
(and invariance of equations). We shall always consider the case when Rq > B, the outer cloak 
boundary. Thus m R > B the solution can be written as ()A.2p . noting that the scattering 
coefficients are equivalent to scattering coefficients for a cavity of radius A. Therein resides 
the reduction in scattering. In a < R < B, the total field is given by Wi + Wg but with an 
argument given by 

R{r) = Co 1 (r + " 2coCi) (50) 

i.e. that corresponding to the hyperelastic deformation described above (see eq. (I37p ) . 

We take 30 terms in the modal sum (1A.2P for the wave field, sufficient for convergence of the 
solution. Figure [2] shows both the total (top) and scattered (bottom) fields corresponding to 
the following problems: scattering from a cavity of radius A with KsA = 27r in an undeformed 
medium (left) and scattering from a cavity with the presence of a hyperelastic cloak (right) 
with undeformed (A) and deformed (a) inner radii defined via a = (3 A where /3 is defined 
above. The outer cloak boundary B is defined by KgB = 4tt (right). This demonstrates 
significantly reduced scattering due to the presence of the hyperelastic cloak as compared 
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with the non-cloaked case. Indeed we are able to quantify this by determining the reduction 
in scattering cross-section, defined in ()A.3P for plane wave incidence. Without the cloak we 
have ^shKs = 5.39 whereas with the cloak -ysnKs = 2.61 resulting in a 51.5% reduction 
in scattering. Figure [3] shows the scattering cross-section ^shKs (left) together with the 
percentage reduction in scattering (right). 
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Figure 4: Scattered fields for the in-plane P/SV problem for an incident field generated by a 
compressional source at Rq = Svr, Go = 0. Left: Compressional (top) and shear (bottom) fields 
corresponding to an undeformed cavity with scaled radius KpA = 27t. Right: Compressional 
(top) and shear (bottom) fields corresponding to a conventional cloak generated via pre-stress 
where the scaled deformed inner radius is KpU = 2tt and initial inner cavity radius is KpA = vr 
so that a = (3 A where /3 ~ 1.863. 

6.2 P/SV wave propagation 

In the P/SV case, in addition to the shear wavenumber Kg we also introduce the compressional 
wavenumber Kp via = Pquj /(Ao + 2/io). We use the undeformed medium solution 

as derived in Appendix IA.2I as a means of determining the solution for the cloak problem. 
This solution is employed in the exterior region together with the same solution but with 
modified argument (due to the hyperelastic deformation) in the cloak region. Thus m. R> B 
the solution can be written as (IA.5P with scattering coefficients An and -B„ given by (TOa]) 
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and ()A.6bp respectively, noting that they are equivalent to scattering coefficients for a cavity 
of radius A and therefore a reduction in scattering is present. Note that here a different effect 
is introduced as compared with the SH case: shear waves are produced as a result of mode 
conversion on the boundary of the cavity. In a < R < B, the total field is given by the 
sum of the scattered and incident fields but with the argument as given in (|50p due to the 
hyperelastic deformation. 

We take 30 terms in the modal sums (jA.Sp . which is sufficient for convergence of the 
solution. Figure m shows the scattered fields corresponding to the P-wave (top) and S-wave 
(bottom) fields associated with = 1/3 and for the following problems: scattering from a 
cavity of radius A with KsA = 27r in an undeformed medium (left) and scattering from a 
cavity with the presence of a hyperelastic cloak (right) with undeformed (A) and deformed 
(a) inner radii defined via a/ A = f3. The outer cloak boundary B is defined by KgB = Att 
(right). Scattering is significantly reduced due to the presence of the hyperelastic cloak as 
compared with the non-cloaked case although it is relatively difficult to see this directly with 
the plots. As with the SH case, let us quantify this by determining the reduction in scattering 
cross-section, defined in (jA.SP for plane wave incidence. Without the cloak ■jpKp = 13.564 
whereas with the cloak "fpKp = 7.258 resulting in a 46.48% reduction in scattering. Figure 
[5] illustrates the scattering cross-section ^pKp (left) together with the percentage reduction 
in scattering (right) compared for three different Poisson ratios: u = 1/3,7/15 and 49/99. 
Note that for very low frequencies there is a huge reduction in scattering, close to 100%. This 
tails off at higher frequencies but still remains at around 50% reduction in scattering which 
is clearly significant. Reduction is larger for smaller Poisson ratios. We also note the rather 
interesting result that the peak in the cross-section actually induces an increase in scattering 
at some values of Kpa as compared with the case without the cloak although this is only for 
a narrow range of such values. This can be associated with the increasing disparity in the P 
and and SV wave numbers as u tends to \, noting that K^/Kl = 2(1 - i/)/(l - 2v). 

7 Conclusions 

The close correspondence between transformation elasticity and small-on-large theory points 
to a method for realizing the former. Specifically, the semilinear strain energy function 
of eq. p9p yields the correct incremental moduli required for transformation of isotropic 
elasticity. The connection between the two theories is that the transformation equals the 
finite deformation. The fact that the pre-stress must be in a state of equilibrium places 
a constraint on the type of transformations allowed. Specifically, they are limited by the 
condition (j24p . or equivalently, trV =constant, which yields stable radially symmetric pre- 
strain |29] . This implies that the actual size of a cylindrical target can be increased in area by 
a factor of 4, its radius by factor of two, without any change to the scattering cross-section. 
The restricted form of the transformation is not surprising considering the fact that the theory 
can simultaneously control more than one wave type, in contrast to acoustics. 

In the two-dimensional problems for which results were provided, it was shown that the 
presence of a conventional cloak generated by the use of pre-stress leads to a significant 
reduction in the scattering cross-section from the cavity, as compared with scattering from 
a cavity without a cloak. This effect is particularly striking at low frequencies and for small 
Poisson ratios. We should note that in general one has to consider stability of nonlinear elastic 
solids in the large deformation regime. While we have not undertaken a full stability analysis, 
we have noted that the deformation for what we have termed hyperelastic cloaking (HC) is 
automatically stable (see ^5.2.3p . Extension of these results will be the subject of subsequent 
study. We also note that manufacturing nonlinear elastic solids with specific strain energy 
functions can be difficult to achieve in practice, although this is certainly no harder than 
generating complex metamaterials which appears to be the current alternative. 

This work sheds some light on transformation methods in other wave problems. In acous- 
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Figure 5: Scattering of P/SV waves from a cylindrical cavity. Left: Scattering cross section 
from the undeformed cavity (solid) with radius Kpa and from a deformed cavity with initial 
scaled radius KpA such that a = (3A where (3 ~ 1.863 (dashed). Right: Percentage reduction 
in scattering cross section 7% due to pre-stress. We have = 1/3 (top), = 7/15 (middle) 
and V = 49/99 (bottom). Note that for the latter case the peak in scattering cross section 
results in a narrow range of values of Kpa where the cloak increases scattering. For other 
values, there is significant reduction in scattering, especially at very low frequencies. 



tics and electromagnetism there is no constraint on the transformation; any one-to-one map- 
ping is permitted. In principle, there is no constraint for transformation elasticity either, 
although the transformed materials are quite difficult if not impossible to obtain, especially 
since they are required to lose the minor symmetry in their corresponding elastic modulus 
tensor. The equivalence of transformation elasticity and small-on-large theory provides a 
unique and potentially realizable solution, although with a limited range of transformations 
allowed. It would be desirable to relax this constraint, which interestingly, does not appear 
for the related problem of SH wave motion in incompressible hyperelastic solids |20j . The 
limit of incompressibility offers a clue to a possible resolution for solids with Poisson's ratio 
close to one-half, and will be the subject of a separate study. 



16 



Appendix 

A Elastic wave scattering from cylindrical cavities 

Brief summaries of the two wave scattering problems are given below. For further details 
see e.g. Eringen and Suhubi [33]. Scattering is considered from a cylindrical cavity of radius 
A, located at the origin of a Cartesian coordinate system X = (Xi, X2, X3), related to 
a cylindrical polar coordinate system via X = {RcosQ, RsinQ, Z). An incident wave is 
generated by a line source of amplitude C (a force per unit length) located at the point 
(i?0) ©o)- We take ©o G [0, 2tt), subtended from the positive X axis. 



A.l SH wave scattering 

In this case the line source is polarized in the Z direction thus creating incident horizontally 
polarized shear (SH) waves which are then scattered from the cavity without mode conversion. 
The total wave field in this domain will therefore be U = (0, 0, W{X, Y)) where W satisfies 

(V^ + K^,)W = CRo~^6{R - Ro)6{Q - Go) (A.l) 

with Kg = pu'^/no and C = Cq/^q. We seek W in the form W = Wi + Wg where 
Wi = (C/4i)Ho(K.5) is the incident field and S = ^/{X - X^^Y + {Y - Yof with Xq = 

i?ocoseo,lo = i^osinGo. We have defined Ho(K,5) = h[)^^(K,5) = Jo(K,5) + iYo(A',5), 
the Hankel function of the first kind, noting that Jq and Yq are Bessel functions of the first 
and second kind respectively, of order zero. Together with the exp(— iwi) time dependence in 
the problem, this ensures an outgoing field from the source. Graf's addition theorem allows 
us to write this field relative to the coordinate system (i?, 0) centred at the origin of the 
cavity [33] and we can use the form appropriate on i? = ^ in order to enforce the traction 
free boundary condition ^QdW/dR = on i? = ^, yielding the scattered field in the form 



VenI)„H„(i^,fl)cos(n(e-eo)) with = C--p-^Bn{KM. (A.2) 



n=0 



where H„ and J„ are respectively Hankel and Bessel functions of the first kind of order 
n. We have also defined = 1, En = "2, n > 1. If we take Rq ^ 00 and Co = 
2i//o\/2vrAJ2^e*(''/^-'^'"f^o), the incident wave of unit amplitude takes the (plane-wave) form 
Wi = exp{iKs{X cos Qinc + ^sinGi,ic)} where @inc = Oo — £ [— vr, vr) is the angle of in- 
cidence subtended from the negative X axis. The scattered wave Ws takes the form ([A^ i 

with Dn — )• Dn"^^ = —i^y{KsA)/i{'^{KsA). The scattering cross-section of the cylindrical 
cavity for plane wave incidence is [35j 

00 

n=0 



A.2 P/SV wave scattering 

In this case the line source at {Rq, ©o) with amplitude Co is a compressional source. Thus the 
incident field consists purely of in-plane compressional waves. Due to mode conversion, the 
scattered field consists of coupled in-plane compressional (P) and vertically polarized shear 
(SV) waves. The total wave field will therefore be U = (U{X,Y),V{X,Y),0) and using the 
Helmholtz decomposition U = V<I> + V x (^k), we deduce that 

V'^<^ + K^^ = CRo~^6{R-Ro)S{e-eo), V^^- + = (A.4) 
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where = uj'^p/{\q + 2;Uo)) K'l = ^^^p/z^o and C = Cq/{\q + 2ij,q). Seek the wave field in 
the form ^ = + ^s, ^ = where $j = (C/4i)Ho(i^'pS) is the incident compressional wave 
with notation defined in Appendix lA.ll We satisfy the traction free (cjj/j = 0, a^Q = 0) 
boundary condition on i? = A, by using Graf's addition theorem, and the scattered field is 

oo oo 

= Y,£nAnHn{KpR)cos{n{Q - Go)), = Y,enBnHrXKsR)s\n{n{Q - Go)). (A.5) 

n=0 n=0 

The scattering coefficients are 

An='- CYiniKpRo) [xi{KpA)Mf{K,A) - Xl{KpA)Ml^{KsA)\ /A„, (A.6a) 
Bn = \ CYiniKpR^) [Xl{KpA)M}^\KpA) - Xl{KpA)Ml\KpA)] /A„, (A.6b) 

where 

Xl{x) = {n^ + n - \{KsAf)Jn{x) - xJn-i{x), 

X^{x) = n{n + 1) J„(x) - nxJ„_i(x), 
M^\x) = -Mf (x) = {n^ + n- \{KsAf)Hn{x) - xHn-i{x), (A.7) 
Mni^) = -Mn{x) = -n(n + l)//„(x) + nxHn^xix), 

A„ = M^^(KpA)Mf(K,A) - Mf(KpA)Mf(K,A). 

If we take i?o ^ oo together with Co = 2i(A + 2^o) ^JI-KK^RQe^'-'^ 1^-^^^*^'^ the incident wave 
of unit amplitude takes the (plane-wave) form = exp{zi^p(X cos 0j„c + ^ sin G^c)} where 
Gj„,c is defined above in Appendix lA.ll The plane wave scattered fields take the form in 
foil with An, Bn A^r?'"\ B^r?""^ where the latter are given defined in (IA.6P under the 
replacement jClin{KpRo) — )• —i"-. The scattering cross-section 7p of the cylindrical cavity 
for plane compressional wave incidence (subscript P indicating this fact) is |35j 
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